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ABSTRACT 

The arching phenomenon is an emergent pattern formed by a 
c-sized crowd of intelligent, goal-oriented, autonomous, het¬ 
erogeneous individuals moving towards a w-wide exit along 
a long IT-wide corridor, where W > w. We collected empir¬ 
ical data from microsimulations to identify the combination 
effects of c and w to the time T of the onset of and the 
size S of the formation of the arch. The arch takes on the 
form of the perimeter of a half ellipse halved along the minor 
axis. We measured the S with respect to the lengths of 
the major M and minor m axes of the ellipse, respectively. 
The mathematical description of the formation of this phe¬ 
nomenon will be an important information in the design 
of walkways to control and easily direct the flow of large 
crowds, especially during panic egress conditions. 

1. INTRODUCTION 

A crowd is a large group of humans, usually numbering by 
the hundreds to tens of thousands, who are in the same phys¬ 
ical environment, such as a hallway, sharing a common pur¬ 
pose, such as exiting through a sole door along the hallway. 
Although the individuals in the crowd are sharing a common 
purpose, they may act differently in a crowd than when they 
are by themselves Ea- The general movement of a crowd is 
affected by a number of different individuals who are moving 
independently from each other, but sharing a common envi¬ 
ronment, and aiming to reach their respective destinations. 
The understanding of crowd movement is important to plan¬ 
ning and improving shared public places, not only to effec¬ 
tively and efficiently facilitate the comfortable movement of 
individuals, but also to guarantee the individual safety, espe¬ 
cially under conditions of danger when quick and orderly 
evacuation of a mass of individuals is desired. 

The patterns of behavior of individuals in a crowd is affected 
directly by the interactions and influences between the indi¬ 
viduals and their environment. Because data from real evac¬ 
uation is hard to measure and conducting experiments on 
humans has ethical questions, the characterization of the 


dynamic aspects of evacuation processes have been confined 
by experts (engineers, sociologists, and scientists alike) only 
through simulation and modeling. Because of this, several 
computer models have been proposed and utilized which 
aided in the understanding of crowd behaviors under normal, 
controlled egress, and panic situations. Of course, these 
models have been separately evaluated either by comparing 
the existing descriptive real-world data usually compiled by 
social scientists with the descriptive behavior of the model 
outputs, or by comparing the quantitative and qualitative 
outputs of the model with that of another model. One 
among the many models of crowd dynamics is the social 
force model (SFM) Ea which describes the motion of an 
individual as affected by both physical and psychological 
forces. SFM has been a popular model among researchers 
and planners because it was successful in describing phe¬ 
nomena observed in the real world, such as the “faster- 
is-slower” in escape panic and the “arching” phenomenon 
observed as a side-effect to “clogging” on exit ways. How¬ 
ever, when using SFM for crowds under low crowd density 
conditions, empirical results show that the simulated indi¬ 
viduals do not behave as expected. In other words, the 
simulated individuals act irrationally rather than rational 
ones, as we expect of humans under normal conditions. For 
example, the simulated individuals go directly and repeat¬ 
edly towards an object that blocks an exit way, instead of 
finding a way to go around the object and exit the hallway. 
This happens because the desired direction of the simulated 
individual is always towards the exit points, even in the pres¬ 
ence of obstacles between its current position and the exit 
points. 

The SFM has been used extensively by many researchers Ei- 
E1 because of its capability to simulate the often observed 
human crowd phenomena such as flocking and bidirectional 
lane formation, among others. However, the simulated 
individuals under SFM do not exhibit some of the common 
human behaviors observed in a large crowd. These behaviors 
are: 


1. Imitation. According to Festinger @|, there is a drive 
in a person towards its own self evaluation. This drive 
forces the person to belong to groups so that he can 
be associated with others. Thus, in a crowd of moving 
people, the individuals tend to move into groups whose 
members they believe would have the same opinion as 
theirs. 


2. Contagion. Using the imitation behavior, people will 
also tend to “adopt” the behavior of others in the same 
group. Thns, if a group member senses that there is a 
discrepancy in a group with respect to behavior, opin¬ 
ions or abilities, that same member will reduce the 
discrepancy by adcmting the most common behavior, 
opinion, or ability [g]. 

Because of these, researchers 0 recently proposed the 
social comparison theory (SCT) that algorithmically mimics 
the above-mentioned behaviors and was introduced as an 
alternative improvement to SFM. The key idea of SCT 
is that humans who lack the objective means to evaluate 
their current state, such as when in confusion, would com¬ 
pare themselves to others that they perceive are similar to 
themselves (i.e., both imitation and contagion behaviors). 

In our previous work [^, we used the SCT to design an arti¬ 
ficial individual whose behavior resembles that of a human 
when walking (or traveling) within a crowd. We included 
in the same model a “field of desire” to allow the simulated 
individual to change its desired direction of movement volun¬ 
tarily, based on its perception of the behavior and location 
of other individuals and obstacles along the hallway. The 
“field of desire” is a 100° field of vision facing the forward 
motion of the simulated individual, which the individual 
uses to obtain spatial information concerning its surround¬ 
ings and to choose the desired forward direction towards the 
exit goal. Our crowd of simulated individuals were able to 
exhibit the arching, clogging, and bursty exit phenomena 
usually observed in real crowd situations. The arch exhib¬ 
ited at the edge of the clogged crowd resembles that of a half 
ellipse, with the major axis parallel with the direction of the 
crowd flow, while the minor axis is coincident with the axis 
of the exit width. 

We present in this paper the result of the study we conducted 
to understand the interactive effects of crowd size c and exit 
width w in the formation of arch in a clogged crowd moving 
towards one exit door along and in the middle of a hallway. 
We will show by empirical means that the onset of arching 
has a negative correlation with the increasing w, while its 
slope tends to grow steeper as the c is increased. We will 
also show that the profile of the arch resembles a half ellipse 
whose major axis runs along the length of the hallway, while 
its minor axis is parallel to the width of the exit door. The 
major axis grows smaller as the exit width is increased, but 
grows bigger with the increase of the crowd size. The minor 
axis, on the other hand, grows bigger either as the crowd 
size is increased, or as the exit width is increased (of course 
limited only by the width of the hallway). 

2. REVIEW OF LITERATURE 

The approaches used by researchers for crowd simulation 
maybe classified into two categories namely, macroscopic 
and microscopic. The macroscopic models are models that 
describe the dynamics of observable quantities of a crowd 
motion, such as density and velocity Examples of such 
models are the ones used by Hughes [U], Huang et al. [^ . 
and Treuille et al. wherein they captured the density and 
velocity of the crowd using partial differential equations, the 
usual mathematical descriptions of the dynamics of fluid. 


Similar to what Ngoho and Pabico did in the past, 
Musse et al. proposed a model based on the trajecto¬ 
ries of motion captured automatically from video sequences. 
While Ngoho and Pabico [2^ captured the real-time trajec¬ 
tory of the velocity field (see for example Figure [T]), Musse 
et al. used an offline method to cluster similar trajecto¬ 
ries to extrapolate the velocity field. 

Borrowing some ideas in gas-kinematic models, Helbing 
and Molnar [l^ introduced the SFM to simulate pedes¬ 
trian flows, wherein a self-driven particle interacts through 
social rules and regulations, moves in its desired speed and 
direction, and attempts not to collide with obstacles, other 
particles, and surrounding barriers. Thus, in order to reach 
its destination faster, pedestrians take detours even if the 
taken route is crowded [T^. However, the choice of detour 
is dependent on the recent memory of what the traffic was 
like the last time they took the route, the figure of which 
was found by Ganem to be polygonal in nature. In 
agreement with the SFM, Weidmann observed that, 
as long as it is not necessary to go faster, such as going 
down a ramp, a pedestrian prefers to walk with his or her 
desired speed, which corresponds to the most comfortable 
walking speed. However, Weidmann further observed 
that pedestrians keep a certain distance from other pedes¬ 
trians and borders. The distance between the pedestrians 
decreases as the density of the crowd increases. The pedes¬ 
trians themselves cause delays and obstructions. Arns @] 
observed that the motion of the crowd is similar to the 
motion of gases and fluid, while fl8|| suggested that it is 
similar to granular flow as well. 

Helbing et al. [l^. in his extension of the SFM, showed 
that many aspects of traffic flow can be reflected by self- 
driven many-particle systems. In this system, he identified 
the various factors that govern the dynamics of the particles 
such as the specification of the desired velocities and direc¬ 
tions of motion, the geometry of the boundary profiles, the 
heterogeneity among the particles, and the degree of fluc¬ 
tuations. One such observable pattern is the formation of 
lanes of uniform walking direction, formed because of the 
self-organization of the pedestrians [l5||. Aside from the 
self-organizing behavior of the crowd, obstacles were also 
observed to both positively and negatively contribute to the 
flow of the traffic. During escape panic of large crowds, sev¬ 
eral behavioral phenomena were observed [IJ: build up of 
pressure, clogging effects at bottlenecks, jamming at room 
widening areas, faster-is-slower effect, inefficient use of alter¬ 
native exits, initiation of panics by counter flows, and impa¬ 
tience. It was observed that the main contributing behaviors 
in these situations is a mixture of individual and grouping 
behavior. 

While macroscopic models describe the dynamics of the col¬ 
lective crowd, microscopic models describe the speed and 
position of each crowd member in a given time. Historically, 
Reynolds [U] is known to have pioneered the simulation of 
crowd in this manner. In his work, he simulated a flock 
of birds where each simulated bird is implemented as an 
independent individual that navigates according to its local 
perception of the dynamic environment. In a similar study 
in the years that follow, Tu and Terzopoulos simulated 
schools of artificial fishes using perception to drive each fish’s 
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Figure 1: An example output frame of the system developed by Ngoho and Pabico [5^. The upper left 
part contains the input image of student pedestrians with trackers (bounding boxes). The upper right part 
contains the trajectory of the pedestrians. The lower part contains the distance vs. time graph of the 
trajectory. Thi s figure is in color in the electronic copy of this paper printed with permission from Ngoho 
and Pabico [2^ and the Philippine Society of Information Technology Educators Foundation, Inc. (PSITE). 



locomotion while being constrained by underwater physics. 
The then novel programming of the fish’s motion and envi¬ 
ronment resnlted into observed real-world school of fishes’ 
behaviors. 

Because of the success of microsimulation in exhibiting 
observed behaviors in animals, researchers in varions fields 
of science (e.g., transportation science, computer gaming, 
artificial intelligence, and industrial robotics) adopted the 
programming method to simulate the motion and envi¬ 
ronment of human pedestrians in a crowd. The main 
purpose is to aid them in understanding some emergent 
crowd behaviors so that they can design better and safer 
walkways, or build realistic gaming elements, environments, 
and scenarios. During the advancement of the method 
itself, several enhancements to the model, as well as alter¬ 
natives, have been suggested. For instance, Kapadia et al. 
[ 2 ^ proposed the use of what they called an “egocentric 
affordance field” to obtain the speed and the direction of 
each simulated pedestrian. On the other hand, Varas et al. 
[ 3 ^ and Perez et al. introduced cellular automata (CA) 
models to direct the motion of the simulated crowd member. 
Both works located their crowd members at cells of a grid 
wherein their respective coordinates are updated at discrete 
time intervals. The implementation of the CA model was 
so simple that Teknomo implemented both the CA 
simulation and visualization using only a desktop electronic 
spreadsheet program. The simplicity of CA prompted other 
researchers [l3 . Im . to further improve it using the prob¬ 
abilistic and statistical functions used in lattice-gas models 
to direct the motion of the simulated member in a cell. The 
lattice-gas model is known to be a special case of CA. 

Although the above models were already successful in sim¬ 
ulating observed phenomena in crowds, the behavior of the 


crowd member as an individual has not been fully inte¬ 
grated into these models. Because of this, Lerner et al. 
proposed an approach for fitting human behaviors, such as 
talking on the phone, combing hair, and moving the head, 
into the simulated individual. In an attempt to simulate 
the collective behavior of humans in groups, Musse and 
Thalmann [2^ used a hierarchical crowd organization with 
autonomous groups to simulate the grouping or flocking 
behavior in pedestrians. To be able to model and generate 
low-level behaviors for the simulated individuals. Golden- 
stein et al. [ 1 ^ presented non-linear dynamic systems to 
model the interaction of the simulated individuals to its 
environment. Such low-level behaviors are real-time target 
tracking and obstacle avoidance in a highly dynamic envi¬ 
ronment. Later, the authors [l^ used a hybrid low-level 
system integrating the non-linear dynamical systems with 
kinetic data structures to simulate a complex crowd envi¬ 
ronment that includes three-dimensional steering, crowding, 
and flocking within moving and non-moving obstacles. 

Recently, established works in the field of social psychology 
have been successfully incorporated into microsimulation 
models of the crowd, resulting into a much more realistic 
simulated collective behavior. Example of works go as far 
as 1895 such as that of Le Bon [^, and the newer ones 
such as that of Allport [J|, Blumer [3|, and Berk [^. In 
these works, the authors discussed that people in crowds act 
similarly to one another, usually in a coordinated fashion 
that when viewed from a distance, it would appear that 
they are being choreographed by a single entity. However, in 
reality, the coordinated motion is achieved with little or no 
verbal communication at all. This somewhat homogeneous 
motion by the crowd is explained by Le Bon as being 
directed by two processes, imitation and contagion that 
we discussed earlier in Section [T] The coordinated motion 








occurs in a “circular reaction” process, a type of interstimu¬ 
lation between two individuals a and b where the response 
of h to the stimuli from a reproduces a stimulation that is 
reflected back to a. The reflected stimulation reinforces o’s 
stimuli in return, thus resulting into a circular reaction 
process [^. Following the adage that “similar people act 
in similar ways,” Allport explained that the somewhat 
homogeneous motion of the crowd is the product of like- 
minded individuals. The likemindedness of the individuals 
in the crowd has been given a game-theoretic perspective 
by Berk [^. Here, he said that the coordinated behavior of 
the crowd is consistent with that of an individual applying 
a minimax strategy. Given an action A being performed 
by the greatest number of individuals in a crowd, an indi¬ 
vidual c will see A as the action with least cost, and thus 
will perform it also. The individual that selects the action 
of the majority has to pay the least cost compared to that 
which follows the minority. 

One of the most recent models that incorporates works in 
social psychology, includin g th at of cognitive sciences, is that 
of Fridman and Kaminka [10] where they used the ideas of 
SCT by Festinger Qj. In their work, they proposed an algo¬ 
rithmic framework for SCT which resulted into simulated 
crowd behavior that is more in tune with observed real-world 
crowd behavior compared to those techniques that do not 
incorporate the cognitive models. Because of this, we have 
decided to start with the SCT, but hybridized it with our 
own algorithm for mapping the perceived trajectory of an 
individual’s goal, limited W the individual’s held of vision 
(i.e., our “held of desire”) [0. Our hybrid model was able to 
exhibit important crowd behaviors: arching, clogging, and 
bursty exit rates. We have also observed a new phenomenon 
we called double arching, which we will discuss in detail in 
another paper and forum. In this paper, we will explore 
the interactive roles of c and w to the formation of the arch 
using the hybridized microsimulation approach in the hope 
of exploiting the outcomes for designing better and safer 
walkways under various crowd egress situations. 

3. METHODOLOGY 

3.1 Programmin g the Individual 

In our approach, each simulated individual was programmed 
with the following general behaviors which were described 
originally by Wooldridge and Jennings [3^ . and then later 
on extended and explained further by Epstein |3, Franklin 
and Graesser [^, Torrens [s^ and Macal and North |^: 

1. Autonomy. The individuals are autonomous units, and 
thus are not directed by any centralized control. They 
are capable of processing information and exchanging 
this information with other individuals in order to 
make independent decisions. They may freely interact 
with others in a few given situations, and this freedom 
does not affect their being autonomous. With this, we 
say that the individuals are active rather than purely 
passive (see discussion of “active” below). 

2. Heterogeneity. We believe that individuals should be 
heterogeneous to allow the microsimulation to develop 
autonomous individuals. 


3. Active-. Being active means that the individuals influ¬ 
ence a simulation independently. The following active 
features can be identified: 

(a) Pro-active or goal-directed 

(b) Reactive or perceptive 

(c) Bounded rationality 

(d) Interactive and communicative 

Aside from the abovementioned general behaviors, we infuse 
the SCT by following the algorithmic approach of Fridman 
and Kaminka m, which we describe here for clarity. The 
SCT follows a subset of axioms as follows: 

1. When individuals can not objectively decide what to 
do next, they compare their current state with that of 
others. Each individual Vi is characterized by a tuple 
with k attributes V — {of, 02 ,... , a^), where each 
attribute a) of individual Vi, VI < j < k, corresponds 
to a point in a fc-dimensional space. Each dimension 
corresponds to an attribute such as the coordinates of 
Vi in two-dimensional space, heading, height, gender, 
etc. 

2. Individuals compare their state to individuals that are 
similar to them. We defined a similarity function 
Sim{Vi,Vj) which measures the similarity between 
the ith and jth individuals (i 7 ^ j). The similarity 
function is given in Equation [T] For an individual Vi, 
we selected the individual with the highest similarity. 

3. If there are differences in their state with the object 
of comparison, the individuals take steps to reduce the 
differences. If the similarity of the selected individual 
is below a given value, then Vi is triggered to do some¬ 
thing to reduce the difference. 

We measured the similarity between Vx and Vy indepen¬ 
dently along each dimension, such that the similarity in 
dimension Oi is a function Sai{af,af) : Oi x Oi [ 0 , 1 ]- 
A value of 1 indicates full similarity, while a 0 means dis¬ 
similarity. For example, one attribute of Vx is its distance 
from Vy, thus a value of 1 means they are positioned in the 
same location, while a 0 means they are far apart. 


k 

Sim{Vx,Vy) = '^Sa,j{aj,a^.) X Wj (1) 

Because our simulated individuals exhibit these behaviors, 
our simulation approach is more akin to modeling humans 
and objects in very realistic ways than other modeling 
approaches, such as those that aggregate mathematical 
equations to explain the dynamics of pedestrian movement 
(i.e., macrosimulations). Thus, our simulation can be used 
with higher confidence to perform what-if scenarios to aid, 
for example, university administrators and decision makers 
with regards to management policies, as well as infrastruc¬ 
ture development plans, for safer learning environments for 
the students and constituents (see for example Figure [2]). 



Figure 2: Snapshots of example what-if scenario simulations on the effect of exit door configurations to student 
egress (students are shown as black circles with pointers to visualize heading) in a very large lecture hall: 
(a) Current exit door configuration which opens outward but hinders student flow; (b) Effect of moving the 
door hinge such that the exit door opens behind the student flow; and (c) Effect of increasing the effective 
width of the exit door by 50% and replacing the swinging door by a sliding door. This figure is in color in the 
electronic copy of this paper printed with permission from Castro and Pabico and the Philippine Society 
of Information Technology Educators Foundation, Inc. (PSITE). 


3.2 The Goal of the Individual 

In a two-dimensional flat world, let G be the goal of Vj (Vj) 
with set coordinates {{xi,yi),{x2,y2),. ■ ■ ,{xn,yn)} repre¬ 
senting a line segment. This set is nsually the location of a 
fixed-width exit door. Since all V aim to reach any one of the 
exit coordinates, a given Vi will face towards the direction 
of the nearest exit coordinates from its cnrrent coordinates. 
The Vi will scan its field of vision for the closest free space 
and move towards that space with its gaiting speed, which 
is currently set at one pace per simulation time step. In the 
future, we wish to vary this gaiting speed depending on the 
height of Pi. A free space is a location in the flat world 
that is not a wall, a Vj where j yf i, or any other object. 
If other individuals are blocking Vi’s space within its held 
of vision, Vi stops at its current coordinates. It is possible 
that the individual might move away from the target exit 
coordinate if the chosen free space was at the edge of its 
held of vision. When that happens, the individual will still 
move to the free space, but after moving, it will redirect its 
heading towards the possibly new nearest exit coordinates. 
When the distance of the individual to the nearest exit coor¬ 
dinates is < 1, it considers itself as already exited and will 
move to the edge of the simulation world. Until this dis¬ 
tance has not been achieved, the individual will just repeat 
its decision making process as described above. 

3.3 Interactive Effects of Crowd Size and Exit 
Width on Arching Profile 

We conducted experiments along a fixed-width hallway to 
determine the interactive effects of c and w to the time T 
of the onset of formation of the arch, as well as the dimen¬ 
sion of the arch’s profile at T. The arch profile is a half 
ellipse, halved across the major axis, and is measured with 
respect to the length of the minor and major axes. We 
experimented on five levels of c (200, 300, 350, 400, and 
450 individuals) and seven levels of w (1, 3, 5, 7, 9, 11, 
and 13), resulting into 35 c-w combinations. We replicated 


each combination three (3) times to take into considera¬ 
tion random variabilities brought about by external factors 
that can not be modeled, such as the implementation of the 
random number generator of our host operating system and 
other systemic variabilities such as the occurrences of hard¬ 
ware interrupts during simulation that might have affected 
our measurements. The exit width is defined as the number 
of individuals in a straight line that can concurrently and 
comfortably yet tightly fit along the axis of the exit width. 
Similarly, the fixed-width hallway, set at 19, is also defined 
as the number of individuals along a straight line that can 
comfortably yet tightly fit in a straight line parallel to the 
corridor width. 

4. RESULTS AND DISCUSSION 

Figure |3] shows an example profile of a crowd that exhibits 
arching near the exit door. We see here that the arch resem¬ 
bles a half ellipse with the ellipse’s major axis running along 
the direction of the crowd flow and perpendicular to the exit 
width, while the minor axis is parallel to the exit width. This 
result is a confirmation of the results of Epstein Q , Franklin 
and Graesser and Helbing et al. [l8||. 

4.1 Arching Time 

Figure |4] shows the response of the onset of arching time T 
to the combination effect of varying c and w. We can see 
here that across different c, there is a recurrent pattern that 
the T decreases as the w is increased, but the amount of 
decrease is not the same between any c. To better quan¬ 
tify this observed decrease in T at increasing w, we con¬ 
ducted a regression analysis between T and w for each value 
of c. Equations [2] through |6] show these regressions, respec¬ 
tively for c = {450,400,350,300,200}. All coefficients are 
significantly different from zero at a confidence level of a = 
0 .01, while all values are significantly different from zero 
at a = 0.05. This means that even though Equations[5]and|6] 
have low values, these values are non zero. Looking at 

























Figure 3: Profile of the crowd before the exit door as it forms an arch. This is a snapshot of the microsim¬ 
ulation when the c = 400 at w = 7. Notice that the major axis runs along the axis of the corridor, while the 
minor axis is parallel with the axis of the exit door. 


the respective slopes of the regression, we see here that as 
the c is decreased, the slope gets to zero. This means that 
at lower c values, the onset of the formation of the arch is 
becoming independent of the w. The w has a significant 
effect on T at higher c. Generally, however, T oc 


T 450 = 

— l.SOw 

-b 

44.72, R? = 

= 0.97 

(2) 

r4oo = 

—2.04w 

-b 

48.25, R? -- 

= 0.92 

(3) 

Paso = 

-0.54w 

-b 

37.89, R^ -- 

= 0.95 

(4) 

Psoo = 

—0.32w 

-b 

35.49, R^ -- 

= 0.66 

(5) 

r2oo = 

—O.lSui 

-b 

34.39, R^ -- 

= 0.34 

(6) 


4.2 Arch Profile 

Figure [5] shows the response of the lengths of the major and 
minor axes at varying values of c and w. Except for c at 450 
and 400 individuals, the length of the major axis decreases 
as the c is decreased, and as the w is increased, while the 
length of the minor axis decreases as the c is decreased, and 
increases as the w is increased. Mathematically, M oc 
while m oc c ■ w. 

We observed here from the figure that when the w = 9, 
the c = 400 plateaued at m = 19, while with M, the slope 
becomes positive. Similar behavior was observed for c = 450 
but the plateauing and the sudden change of direction of 
the slope for m happened at w = 7. At these c values, the 
individuals have already reached the edge of the corridor 
that m cannot expand anymore. Since m is constrained 
to 19, the arch must accommodate more individuals, thus 
translating into the increase in M. However, had the hallway 
width been way more than 19, the same pattern would have 
been observed for all c. We confirmed this assertion when we 
rerun the simulation for c = {400, 450} and w = {9, 11,13} 
at the same number of replicates. When we plotted the 
mean of each c x w combination across replicates, they fall 
reasonably well within the regression line that we obtained 
earlier, of course within some acceptable variance (i.e., the 
respective means’ ±lx standard deviations intersect with 
the regression line, see Figure [6]). 


5. SUMMARY AND CONCLUSION 

In this paper, we present the behavior for a simulated indi¬ 
vidual used in microsimulation of pedestrian crowd. The 
behavior is based on the SCT and “field of desire” hybrid 
technique in our attempt to provide cognitive behaviors 
to the simulated individual so that the crowd can exhibit 
crowd patterns that are observed in real-life, such as that 
of arching. With this programming, we experimented on 
the combination effect of c and w on T, and the arch pro¬ 
file described by M and m. Here, we found out that the 
T oc the M oc —, and the m oc c x w. These results 
suggest that, intuitively, the w should be wide enough at 
higher c so that T will happen at the least possible time, 
and at the same time M will be shorter, and m will be 
longer (but still M > m). We believe, though, that when 
we find the values for c and w where M = m, then this is 
the combination when T will be the fastest. Thus, we will 
investigate this hypothesis in the future. 
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Figure 4: T as affected by c and w combination: (a) c = 450; (b) c = 400; (c) c = 350; (d) c = 300; and (e) c = 200. 
The blues squares are means of T at each c—w combination, averaged across all replicates. The blue lines are 
the respective regression lines. 



Figure 5: (a) M and (b) m as affected by c and w combination. Colored shapes are means while lines are 
regression fit at each c. Orange/triangle represents c = 450, blue/square c = 400, green/inverted triangle 
c = 350, red/triangle pointing right c = 300, and purple/diamond c = 200. 







Figure 6: (a) M and (b) m as affected by c and w combination when W = 35. 
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